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We propose a simple method to measure synchronization and time delay patterns between signals. 
It is based on the relative timings of events in the time series, defined e.g. as local maxima. The 
degree of synchronization is obtained from the number of quasi-simultaneous appearances of events, 
and the delay is calculated from the precedence of events in one signal with respect to the other. 
Moreover, we can easily visualize the time evolution of the delay and synchronization level with an 
excellent resolution. 

We apply the algorithm to short rat EEG signals, some of them containing spikes. We also apply 
it to an intracranial human EEG recording containing an epileptic seizure, and we propose that the 
method might be useful for the detection of foci and for seizure prediction. It can be easily extended 
to other types of data and it is very simple and fast, thus being suitable for on-line implementations. 

PACS numbers: 05.45.Tp; 05.45.Xt; 87.90.+y; 87.19.Nn 



I. INTRODUCTION 

In recent years, several measures of synchronization 
have been proposed and applied successfully to different 
types of data. Among these studies we can distinguish 
two main approaches: 1) One based on similarities of 
trajectories in phase space (constructed e.g. by time- 
delay embedding) § ||, @, §; 2) One that measures 
phase differences between the signals, where the phases 
are defined cither from a Hilbert Jo, [?], |j| or from a wavelet 
transform H [To) (as shown in ||, these two apparently 
different phases are indeed closely related). 

These new methods compete in popularity with stan- 
dard measures such as the cross-correlation, the coher- 
ence function, mutual information, and also with simple 
visual inspection of the recordings. Cross-correlation and 
coherence are clearly the measures most used so far. In 
contrast to them, all new measures are non-linear in the 
sense that they depend also on properties beyond second 
moments. In addition, some of them have the advan- 
tage of being asymmetric, eventually being able to show 
driver-response relationships || . 

Among others, synchronization measures have been 
used for the study of electroencephalogram (EEG) sig- 
nals. Applications include prediction and localization of 
epileptic activity ]|, ||, |) , phase- locking between different 
recording sites upon visual stimulation [T^] , resonance 
between EEG and muscle activity in Parkinson patients 
, desynchronization upon lesions in the thalamic retic- 
ular nucleus in rats ]Tl| ] , synchronization in motoneurons 
within the spinal cord Q, etc. 
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In the present paper we present a very simple algorithm 
that can be used for any time series in which we can define 
events. These can be spikes in single-neuron recordings, 
epileptiform spikes in EEGs, heart beats, stock market 
crashes, etc. In principle, when dealing with signals of 
different character, the events could be defined differently 
in each time series, since their common cause might man- 
ifest itself differently in each series. This event synchro- 
nization (ES) does not require the notion of phase. It 
cannot distinguish between different forms of m : n lock- 
ings |(| 0], but it can tell which of the two time series 
leads the other. And, above all, it is very simple con- 
ceptually and easy to implement. Due to that, it can be 
used on-line and can show rapid changes of synchroniza- 
tion patterns. 



II. EVENT SYNCHRONIZATION AND DELAY 
ASYMMETRY 

Given two simultaneously measured discrete univariate 
time series x n and y n , n = 1, . . . , N, we first define suit- 
able events and event times if and (i = 1, . . . , m x ;j = 
1, . . . ,m y ). In the signals to be analyzed in this paper, 
these events will be simply local maxima, subject to some 
further conditions. If the signals are synchronized, many 
events will appear more or less simultaneously. Essen- 
tially, we count the fraction of event pairs matching in 
time, and we count how often each time series leads in 
these matches. Similar concepts were used in Jl2[ . 

Let us first assume that there is a well defined char- 
acteristic event rate in each time series. Counter exam- 
ples include strong chirps and onsets of epileptic seizures 
where event rates change rapidly. Such cases will be 
treated below. Allowing a time lag ±r between two 'syn- 



2 



chronous' events (which should be smaller than half the 
minimum inter-event distance, to avoid double counting), 
let us denote by c T (x\y) the number of times an event ap- 
pears in x shortly after it appears in y, i.e: 
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and analogously for c T {y\x). Next, we define the sym- 
metrical and anti-symmetrical combinations 



Qr = 



c T (y\x) + c T (x\y) 



c T {y\x) - c T (x\y) 



which measure the synchronization of the events and 
their delay behavior, respectively. They are normalized 
to < Q T < 1 and —1 < q T < 1. We have Q T = 1 if and 
only if the events of the signals are fully synchronized. 
In addition, if the events in x always precede those in y, 
then q T = 1. 

In cases where we want to avoid a global time scale r 
since event rates change during the recording, we use a lo- 
cal definition for each event pair (ij). More precisely, 
we define 

r« =min{tf +1 -tf,tf-tf_i,*? + i-*?,*?-*"-i} / 2 ' ( 4 ) 

We then define Jy as in Eq.(|J) with r replaced by r^ , 
and c(x\y) as in Eq.(0) with JT. replaced by J i3 . The 
factor 1/2 in the definition of Ty avoids double counting 
if, e.g., two events in x are close to the same event in 
y. Of course, one could also make other choices, e.g. 
by taking Ty smaller than in Eq.(^) or by using — 
min{T, Tij}. As in the definition of events, an optimal 
choice of r depends on the problem. In the following we 
shall suppress the dependence on r, understanding that 
all formulas apply for both variants. 

To obtain time resolved variants of Q and q we simply 
modify eq.(nl) to 



(5) 



with n = 1, . . . , N and 8 the step function (i.e. 0(x) = 
for x < and Q(x) — 1 for x > 0). Similarly, c n (y\x) 
is obtained by exchanging x and y. Then, we define 
the time resolved anti-symmetric combination q(n) = 
Cn(y\x) — c n {x\y) which can be seen as a random walk 
that takes one step up every time an event in x pre- 
cedes one in y and one step down if vice versa. If an 
event occurs simultaneously in both signals or if it ap- 
pears only in one of them, the random walker does not 
move. Exchanging x and y just reverses the walk. For 



non-synchronized signals, we expect to obtain a random 
walk with the typical diffusion behavior. With delayed 
synchronization we will have a bias going up (down) if x 
precedes (follows) y. We should remark that such a bias 
clearly shows the presence of a time delay of one signal 
with respect to the other, but does not necessarily prove 
a driver-response relationship, although it might suggest 
it. In fact, internal delay loops of one of the systems 
can fool the interpretation. Also, the two signals might 
be driven by a common hidden source and the bias just 
indicates different delays. 

The time course of the strength of ES can be obtained 
from Q(n) = c n (y\x) + c n (x\y). If an event is found both 
in x and y within the window r (resp. ry), Q(n) increases 
one step, otherwise it does not change. Of course, Q(n) 
will also not change if there are no new events at all. 
The synchronization level at time n, averaged over the 
last An time steps, is thus obtained as 



Yin) = 



Q(n) - Q(n - An) 



(6) 



where A.n x and An y are the numbers of events in the 
interval [n — An, n]. Similarly, we can also define instan- 
taneous delay asymmetries q'(n). 



III. APPLICATIONS 

Let us now apply these concepts to two sets of intracra- 
nial EEG recordings, one from rats and the other from 
an epileptic patient. 



A. Rat EEGs 

The five pairs of rat EEG signals were obtained from 
electrodes placed on the left and right frontal cortex of 
male adult WAG/Rij rats (a genetic animal model of hu- 
man epilepsy) |l| . They were referenced to an electrode 
placed in the cerebellum, filtered between 1-100 Hz and 
digitized at 200 Hz. In Fig. |l| we show these signals p4[ . 
The first pair (example A in Fig. [l]) is a normal EEG, 
all others contain spike discharges (not to be confused 
with spikes in single neuron recordings) which are the 
landmark of epileptic activity. They arise from abnormal 
synchronization in an epileptic brain even when there are 
no seizures. A localized appearance of spikes can indeed 
delimit a zone with abnormal activity (though this will 
not necessarily be the epileptic focus). Furthermore, time 
delays between them can identify the electrode closest to 
the epileptic focus, especially at the onset of seizures. 

Several measures of synchronization were recently ap- 
plied to the first three cases of Fig. [l] ||. Since spike 
trains lasted usually about 5 seconds, the challenge was 
to try the different measures in these short epochs. Sur- 
prisingly, nearly all the measures gave qualitatively sim- 
ilar results hard to be guessed beforehand. These exam- 
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FIG. 1: Five pairs of rat EEG signals from right and left 
cortical intracranial electrodes. For a better visualization, 
left signals are plotted with an offset. 



pies and two additional cases (D and E), also containing 
spikes, will be further analyzed in this paper. 

For the example A it is difficult, due to its random- 
like appearance, to visually estimate its level of synchro- 
nization and any delay of one electrode with respect to 
the other. However, we can already observe some pat- 
terns appearing nearly simultaneously in both the left 
and right channels, thus showing some degree of interde- 
pendence. The spike-wave trains in the other examples in 
principle suggests a high level of synchronization. How- 
ever, as already shown in ||, the spikes of example C 
appear with a varying time lag between right and left 
channels and are therefore much less synchronized than 
those in B. This is of course not easily seen by visual 
inspection of Fig. [I], but will be clear from the following 
analysis. 

Events were defined as local maxima fulfilling the fol- 
lowing additional conditions: 



1. x(U) > x(t i+ k), for k 

2. x{U) > x{t i±K ) + h 



-K + l, 



,,0, 



K-l 



and the same for y. We took K = 3 and h = 0.1. Other 
choices gave very similar results. 



Example 


Qr = 2 


<?T = 2 


q7=2 


q7=2 


A 


0.57 


0.15 


0.24 


-0.01 


B 


0.80 


-0.29 


0.29 


0.01 


C 


0.48 


-0.20 


0.13 


-0.01 


D 


0.93 


-0.59 


0.41 


0.04 


E 


0.90 


-0.13 


0.46 


0.03 



TABLE I: Time averaged event synchronization Q and delay 
q for the five examples of Fig. |l|. Positive values of q indicate 
that events in the left side lags behind the right one. Surro- 
gate values of synchronization were obtained by shifting left 
channel signals by 500 data points. 



Since the rate of events is more or less constant, we 
used a fixed r. The choice r = 2 gave a good discrimina- 
tion between the five cases. All results shown below were 
compared to those obtained with surrogate pairs which 
were defined by shifting the left channel signals 500 data 
points (2.5 sec) to the right, with periodic boundary con- 
ditions. Our test hypothesis is that without changing the 
individual properties of each signal, after a large enough 
shifting synchronization should reach a background 'zero' 
level. The usefulness of such surrogates was discussed in 
m more detail in [§. 

For the five EEG signals of Fig. [I], we show the values 
of Q T =2 and q T =2 in Table B, both for the original sig- 
nals and the 'time-shifted' surrogates. We observe that 
synchronization levels rank D>E>B>A>C. This is 
in agreement with the analysis of examples A, B and C 
done in p] with several other measures of synchroniza- 
tion. Note that even example A is ranked consistently 
with the other measures, although it does not contain 
obvious events such as the spikes of the other examples. 

All synchronization values are clearly higher than those 
of their corresponding surrogates (surrogates constructed 
with other delay values gave similar results). These sur- 
rogate values vary a lot for the different examples. This 
stresses the importance of keeping the individual proper- 
ties of the signals when constructing surrogates. Except 
for example A, the values of q show that the signals from 
the right hemisphere lag behind the left ones (y). A closer 
visual inspection of Fig. 1 at higher resolution shows that 
this lag is usually 1 data point. The reason of this sys- 
tematic lag is unclear (it could be an artifact of the data 
acquisition or a real physiological effect) and it is beyond 
the scope of this paper. 

In Fig. |^ we show the time evolution of synchroniza- 
tion Q'(n) for the five examples, calculated with a win- 
dow of An = 100 data points. For most of the time, they 
are higher than the values calculated from time-shifted 
surrogates (the light blue horizontal lines indicate time 
averages ±1 standard deviation). In examples A, B and 
C we see abrupt changes of synchronization with time 
which seem statistically significant. In retrospect they 
can also be seen in Fig. [I] on closer inspection, but they 
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FIG. 2: Time resolved event synchronization Q' T=2 {n) for the 
examples of Fig. |l| Blue horizontal lines correspond to the 
time averages ±lcr of the surrogate. 



are much less obvious there and could easily be missed. 
Compared to the first three, examples D and E are more 
stable in time. Finally, the time resolved ES shows a 
better resolution than all synchronization measures con- 
sidered in P]. 

Figure pTshows the time resolved asymmetry between 
the right and the left channels (upper plot) and the re- 
sults from surrogates (lower plot). In all five cases, the 
bias is in agreement with the q values shown in Table |. 
The bias in example D is not only the strongest but also 
the most constant, confirming that D shows the most ro- 
bust and stationary ES (compare Fig. ||). For the other 
examples we see regular changes with time. This is of 
course very difficult to see in the original recordings, 
and it was also not seen with any of the synchronization 
measures studied in 15|. As expected, for the surrogates 
we obtain random walks with small and erratic displace- 
ments. 



B. Human EEG 

As a second example we analyzed an intracranial EEG 
recording from an epileptic patient containing 12 min. of 
pre-seizure and seizure EEG. Data were recorded from 
2 needle shaped depth electrodes with 10 contacts each. 



FIG. 3: Time delays between the right and left channels (up- 
per plot) and for the surrogates (lower plot). Up (down) shifts 
mean precedence of the right (left) channel. 




FIG. 4: Time averaged event synchronization between the 
contacts on the left and right depth electrodes (TL1-10 and 
TR1-10, respectively). 



They were symmetrically placed in the left (contacts TL1 
to TL10) and right (contacts TR1 to TRIO) temporal 
lobes, in the entorhinal cortex and hippocampal forma- 
tion. The EEG was sampled at 173 Hz and band pass 
filtered between 0.53-40 Hz. For further details on the 
data we refer to 0. As in the previous example, event 
times were defined as local maxima, but using K = 10 
and ft. = 50 (this large K was needed because the data are 
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FIG. 5: Time resolved event synchronization (upper plot) and 
delay asymmetries (lower plot) between an channel near the 
epileptic focus (TL7) and the remaining channels on the same 
side. The red bar shows the duration of the epileptic seizure. 



FIG. 6: Delay patterns between an contact in the non-focal 
side (TR7) against the other contacts in the non-focal side 
(upper plot) and against the contacts in the focal side (lower 
plot). No anomalous behavior is seen during the seizure (red 
bar). Notice the different scales in the two plots. 



more noisy than the rat data, and smaller values would 
have led to many spurious events). Due to the varying 
event rate, we used a variable- t approach. For the time 
resolved event synchronization Q'(n) we took a window 
An = 1730. 

Figure ^ shows the time-averaged event synchroniza- 
tion values between all channels. A detailed analysis of 
synchronization patterns for similar recordings has al- 
ready been described by Arnhold et al. Q using a ro- 
bust measure of non-linear synchronization. Here, we 
just summarize the main results which are in perfect 
agreement with those in Ref . || . We first note that syn- 
chronization between left and right electrodes is relatively 
low and that the right contacts form two clusters: TR1-3 
and TR4-10. This is just due to the fact that the first 
3 contacts were located in the entorhinal cortex and the 
remaining ones in the hippocampus ||. Moreover, for 
the right side we observe a gradual decrease of synchro- 
nization with increasing distance between contacts. The 
synchronization pattern for the left channels is different. 
There, the entorhinal cortex/hippocampus separation is 
overshadowed by the epileptic activity leading to a higher 
overall synchronization level. 

A visual analysis of the seizure onset revealed that con- 



tacts TL7 and TL8 showed the first signs of seizure ac- 
tivity. Figure ^| shows the time resolved synchronization 
Q'(n) and delays q(n) between TL7 and the remaining 
left side channels. As expected, synchronization is largest 
between TL7 and its neighbors TL8 and TL6. It is not 
homogeneous in time and we have several short drops 
before seizure starts. Moreover, starting at seizure on- 
set and during the whole seizure, synchronization of TL7 
with TL8 and TL9 is high, while synchronization with 
TL6 and all others is decreased. The lower panel shows 
that all left channels lag behind channel TL7. There is 
just one exception: During the first part of the seizure, 
channel TL7 falls back and channel TL8 leads for about 
half a minute (indeed, the lead of TL7 is weakened al- 
ready some 3 minutes before the seizure). After this, 
TL7 takes up its lead even more vigorously than before. 
This might indicate that the source of epileptic activ- 
ity moves. Whether these features are common to many 
epileptic seizures and whether they can have clinical sig- 
nificance for e.g. seizure anticipation or focus localization 
requires further study with a larger database. 

In Fig. H we show the delays of the contralateral chan- 
nel (TR7) with respect to the other right channels (up- 
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per plot) and to the left channels (lower plot). Channels 
TR4-6 strongly and steadily follow channel TR7, which 
itself follows channels TR8 and TRIO. This might reflect 
the source of 'normal' synchronized activity. A detailed 
analysis is outside the scope of this paper and will be fur- 
ther addressed elsewhere. As seen from the lower panel, 
synchronization between both hemispheres is weak and q 
shows unbiased random walks. The complete absence of 
any deviant behavior during the seizure reflects the fact 
that the seizure does not spread to the contralateral side. 

IV. CONCLUSION 

In conclusion we presented a new approach to mea- 
sure synchronization and time delays that is based on 
the relative timings of events (in this study defined as 
local maxima). This also gives an easy visualization of 
time-resolved synchronization and delay patterns. The 
method is appealing due to its simplicity, straightfor- 
ward implementation and speed. These features make 



very easy its on-line implementation. In the particular 
case of EEGs, the proposed approach is promising for 
the study of recordings of epileptic patients, where syn- 
chronization is important and the analysis of time delay 
patterns could be useful for the localization of the epilep- 
tic focus and the prediction of seizure onset. Also, the 
method should be well suited for single-neuron record- 
ings, where the fast dynamics of spikes makes difficult 
the analysis with other measures. In this paper we fo- 
cussed on application to EEG signals, but the method 
can be easily applied to other types of data just by ad- 
justing the definition of events. 
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